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o 

We study the entropy-corrected version of the new agegraphic dark energy (NADE) model 
and dark matter in a spatially non-flat Universe and in the framework of Hofava-Lifshitz 
cosmology. For the two cases containing noninteracting and interacting entropy-corrected 
I/") ■ NADE (ECNADE) models, we derive the exact differential equation that determines the 

evolution of the ECNADE density parameter. Also the deceleration parameter is obtained. 
Furthermore, using a parametrization of the equation of state parameter of the ECNADE 
model as uj\{z) = cjo + wiz, we obtain both ujq and u%. We find that in the presence of 
interaction, the equation of state parameter coo of this model can cross the phantom divide 
line which is compatible with the observation. 
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i. introduction 

Astronomical observations indicate that our Universe is in a phase of accelerated expansion 
(Perlmutter et al. 1999; Bennett et al. 2003; Tegmark et al. 2004; Allen et al. 2004). One 
explanation for the cosmic acceleration is the dark energy (DE), an exotic energy with negative 
pressure. The dynamical nature of DE, at least in an effective level, can originate from various 
fields, although a complete description requires a deeper understanding of the underlying theory 
of quantum gravity. Nevertheless, physicists can still make some attempts to probe the nature of 
DE according to some basic quantum gravitational principles. Two examples of such a paradigm 
are the holographic DE (HDE) and the agegraphic DE (ADE) models which have originated from 
quantum gravity and possess some of its significant features. The former, that arose a lot of 
enthusiasm recently (Cohen et al. 1999; Hsu 2004; Li 2004; Huang and Li 2004; Jamil et al. 2009a, 
2009b; Jamil and Farooq 2010; Setare and Jamil 2010; Wang et al. 2005a, 2005b, 2008; Sheykhi 
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2010a; Karami 2010a, 2010b; Karami and Abdolmaleki 2010a, 2010b; Karami and Fehri 2010a, 
2010b), is motivated from the holographic hypothesis ('t Hooft 1993; Susskind 1995) and has been 
tested and constrained by various astronomical observations (Feng et al. 2005; Zhang and Wu 2005, 
2007). The latter originated from the uncertainty relation of quantum mechanics together with 
the gravitational effect of general relativity (GR). The ADE model assumes that the observed DE 
effect comes from spacetime and matter field fluctuations in the Universe (Sasakura 1999; Cai 2007; 
Wei and Cai 2008a, 2009). Following the line of quantum fluctuations of spacetime, Karolyhazy 
(Karolyhazy 1966; Karolyhazy et al. 1982, 1986) proposed that the distance t in Minkowski 
spacetime cannot be known to a better accuracy than 5t = where e is a dimensionless 

constant of order unity and tp is the reduced Planck time. Based on the Karolyhazy relation, 
Maziashvili proposed that the energy density of metric fluctuations of Minkowski spacetime is 
given by (Maziashvili 2007a, 2007b) 

where Mp is the reduced Planck mass Since in the original ADE model the age of the 

Universe is chosen as the length measure, instead of the horizon distance, the causality problem 
in the HDE is avoided (Cai 2007). The original ADE model had some difficulties. In particular, 
it cannot justify the matter-dominated era (Cai 2007). This motivated Wei and Cai (2008a) to 
propose the new ADE model, while the time scale is chosen to be the conformal time instead of 
the age of the Universe. The NADE density is given by (Wei and Cai 2008a) 

PA = 2 

where 3n 2 is the numerical factor and rj is the conformal time and defined as 

J a J Ha z 

The ADE models have been examined and constrained by various astronomical observations (Kim 
et al. 2008a, 2008b; Wu et al. 2008; Zhang et al. 2008; Neupane 2009; Sheykhi 2009a, 2009b, 
2010b; Karami et al. 2010; Karami and Abdolmaleki 2011; Wei and Cai 2008b). Besides, in the 
loop quantum gravity, the entropy-area relation can be modified due to the thermal equilibrium 
fluctuations and quantum fluctuations (Banerjee et al. 2008; Banerjee and Majhi 2008a, 2008b, 
2009; Banerjee and Modak 2009; Majhi 2009, 2010; Modak 2009). The corrected entropy takes the 
form 

*=^ + «ln^ + /3, (4) 
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where a and (3 are dimensionless constants of order unity. Taking the corrected entropy-area 
relation (j4| into account, the NADE density will be modified as well. Motivated by the corrected 
entropy-area relation @, the ECNADE density is given by (Wei 2009; Karami and Sorouri 20f0; 
Karami et al. 2011) 

PA = 2^ + - In Miry 2 + tL 5 

where a and f3 are dimensionless constants of order unity. 

Recently a power-counting renormalizable UV complete theory of gravity was proposed by 
Hofava (Hofava 2009a, 2009b). This theory is not Lorentz invariant (except in the infrared limit), 
it is power-counting renormalizable and obeys anisotropic scaling or Lifshitz scaling. Than the 
time coordinate and the 3 spatial coordinates have to be treated separately, the theory is non- 
relativistic, the speed of light diverges in the ultraviolet limit, and test particles do not follow 
geodesies. In consequence causal structures are different from that in General Relativity (Greenwald 
2011) and entropy cannot be defined. Quantum gravity models based on an anisotropic scaling 
of the space and time dimensions have recently attracted significant attention (Visser 2009; Pal 

2009) . In particular, Hofava-Lifshitz point gravity might not have desirable features, but in its 
original incarnation one is forced to accept a non-zero cosmological constant of the wrong sign to 
be compatible with observations (Nastase 2009). There are four different versions of this theory: 
with (or without) projectability condition and with (or without) detailed balance. At a first look 
it seems that this non-relativistic model for quantum gravity has a well defined IR limit and it 
reduces to GR. But as it was first indicated by Mukohyama (Mukohyama 2009, 2010; Saridakis 

2010) , Hofava-Lifshitz theory mimics GR plus dark matter (DM). This theory has a scale invariant 
power spectrum which describes inflation. For reviews on the scenario where the cosmological 
evolution is governed by Hofava-Lifshitz gravity see (Mukohyama 2009, 2010; Saridakis 2010; 
Wang and Wu 2009; Majhi and Samanta 2010; Wang 2011; Khatua et al. 2011). 

Here, our aim is to investigate the ECNADE model in Hofava-Lifshitz cosmology. To do this 
in section 2, we first review the scenario of Hofava-Lifshitz gravity. In section 3, we consider 
a spatially non-fiat Friedmann-Robertson- Walker (FRW) Universe containing ECNADE and DM 
in the framework of Hofava-Lifshitz gravity. We obtain the evolution of dimensionless energy 
density, deceleration parameter and equation of state of ECNADE with interaction/non- interaction. 
Section 4 is devoted to conclusions. 
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II. HORAVA-LIFSHITZ GRAVITY 

Under the project ability condition, the full metric in the (3+1) dimensional Arnowitt-Deser- 
Misner formalism is written as (Calcagni 2009; Kiritsis and Kofinas 2009) 

ds 2 = -N 2 dt 2 + gijidx 1 + N i dt)(dx j + N j dt), (6) 

where the dynamical variables N, TVj and gij are the lapse function, shift vector and 3-dimensional 
metric, respectively. Note that the projectability condition restricts the lapse function N to be 
space-independent, while the shift vector N l and the 3-dimensional metric g^ still depend on both 
time and space. 

Under the detailed balance condition, the gravitational action of Hofava-Lifshitz (HL) gravity 
is given by (Hofava 2009a, 2009b) 

S g = J dtd z x^N ^(KijK ij - \K 2 ) 

t\P ■ • kP fji 

«V /1-4A 
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8(3A- 1) 



-(^R 2 + AR-3A 2 )}, (7) 



where 

^—^(^■-Vi^-VjJVi), (8) 

is the extrinsic curvature and 

C l ' J = e —V k (Ri-- 4 R5j), (9) 

is the Cotton tensor. Also e yfe is the totally antisymmetric unit tensor, A is a dimensionless constant 
and A is a positive constant which as usual is related to the cosmological constant in the IR limit. 
The variables k, w and \x are constants with mass dimensions —1,0 and 1, respectively. Note that 
the detailed balance condition restricts the form of a general potential in a (d + l)-dimensional 
Lorentz action to a specific form that can be expressed in terms of a d-dimensional action of 
a relativistic theory with Euclidean signature, whereby the number of independent couplings is 
considerably limited. 

In particular, Eq.(7) suffers several problems, including instability, inconsistency and strong 
coupling problems, for detail, see (Mukohyama,2010). To overcome these problems, one way is to 
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provoke the Vainshtein mechanism, as first done by Mukohyama for spherical spacetimes (Muko- 
hyama,2010) and by Wang and Wu in the cosmological setting (Wang & Wu,2011). Such con- 
siderations were further carried out by using the so-called gradient expansion method (Izumi & 
Mukohyama,2011;Gumrukcuoglu et al,2011). Another very promising approach is to introduce an 
extra U(l) symmetry, as first done by (Horava &; Melby-Thompson,2010) with A = 1, and later 
generalized to the case with any A by (da Silva,2011). These studies were further generalized to 
the case without the projectability condition (Zhu et al, 2011a; Zhu et al,2011b). In both cases 
(with/without the projectability condition), due to the U{\) symmetry, the spin-0 gravitons are 
eliminated, and all the problems related to them, such as instability, inconsistency and strong 
coupling problems, are resolved. 

To include the matter component we add a cosmological stress-energy tensor to the gravitational 
field equations, by imposing a condition to recover the general relativistic formalism in the low- 
energy limit (Carloni et al. 2009; Leon and Saridakis 2009; Sotiriou et al. 2009; Chaichian et 
al. 2010; Dutta and Saridakis 2010). For this, the energy density p m and pressure p m satisfy the 
energy conservation equation 

p m + 3H(p m + Pm ) = 0. (10) 
In the cosmological context, we use a FRW metric as 

N = l, g lj= a 2 (t) 7lj , N l = 0, (11) 

with 

dr 2 

ltj dx l dxi = - —= + r 2 dn 2 2 , (12) 

1 — kr A 

where k = —1, 0, +1 refer to spatially open, flat, and closed Universe respectively. 

On varying the action with respect to the metric components N and gij, one can obtain the 
modified Friedmann equations in the framework of HL gravity as 

2 - 2 r 3K 2^2 k 2 



H = 77^ TTPm + 



6(3A-l) r 6(3A- 1) L8(3A- l)a 



+ 



3kVA 2 



8(3A-1)J 8(3A-l) 2 a 2! 



(13) 



H + \ R2 = -M&Tf- 

^2 



K 2 p 2 k 2 3K 2 p 2 A 2 



4(3A- 1) L8(3A- l)a 4 8(3A - 1) 
16(3A-l) 2 a 2 ' 



(14) 



6 



where H = ^ is the Hubble parameter. 

Noticing the form of the above Friedmann equations, we can define for DE 

PA = 8(3A - l)a 4 + 8(3A - 1) ' { ' 

K 2 n 2 k 2 3k 2 h 2 A 2 
PA = 8(3A - l)a 4 " 8(3A-1)' ( ' 

The first term on the right hand side proportional to a -4 is effectively the "dark radiation term" , 
present in HL cosmology (Calcagni 2009; Kiritsis and Kofinas 2009), while the second term is 
referred as an explicit cosmological constant. Hence, in expressions f|15j) and (|16p we defined the 
energy density and pressure for the effective DE, which include the aforementioned contributions. 
Finally, note that using (|15p and (|16p it is easy to show that these satisfy the following expression 

p A + 3H(p A +p A ) = 0. (17) 

Finally in order for these expressions to match with the standard Friedmann equations (c = 1) we 
set (Calcagni 2009; Kiritsis and Kofinas 2009) 

K 2 

° c = 16tt(3A-1)' (18) 
KVA 1, (19) 



8(3A- l) 2 

where G c is the "cosmological" Newton's constant. Note that in gravitational theories with the 
violation of Lorentz invariance (like HL gravity) the "gravitational" Newton's constant G g , which 
is present in the gravitational action, differs with the "cosmological" Newton's constant G c , which 
is present in Friedmann equations, unless Lorentz invariance is restored (Carroll and Lim 2004). 
For the sake of completeness we write 

<h = (20) 

Note that in the IR limit (A = 1), where Lorentz invariance is restored, G c and G g are the same. 
Further we can rewrite the modified Friedmann equations (|13[) and ([14p in the usual form as 

TT 9 k 8ttG c , . . 

H 2 + - = —±(p m + PA ), (21) 

(X o 



k +2 R2 + 2^ = - 47rG c(Pm +PA). (22) 
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III. ECNADE IN HL COSMOLOGY 

Here we would like to investigate the ECNADE in HL theory. To do this we consider a spatially 
non-flat FRW Universe containing the ECNADE and DM. Let us define the dimensionless energy 
densities as 

Q Pm _ 8irG c 

— — „ rro Pmi 

Per 3H 2 

PA 87rG c 
Per 3H 2 

k 



n, 



then, the first Friedmann equation (|13|) yields 



i 2 H 2 ' 



(23) 



i — n k = q a + j] m . 



(24) 



From definition p\ = ^ G Oa, we get 



n 



H 2 2^n, 



where 



G c 8irG c 



^ n G g 3n 2 rj 2 



aln(M P 7] 2 ) + [3 



(25) 



(26) 



and M p 2 = 8ttG z . Taking time derivative of Eq. ([5]) and using fj = 1/a, one can get 

2 r Q^2 A/T2 



PA 



at] 



3n z MA a 
2pa + ^ + — 



(27) 



Taking time derivative of Qa = ^jfrPh and using Eqs. (f25j) . ([27|) . rj = 1/a and A = ti^/H, one 
can obtain the equation of motion for £l\ as 



/7 + ^^V /2 (2 7 „-f 



-2^ A n , , 

.ii^ na7„ V 7,, 

8irG c aH 2 Q A 



(28) 



3n 4 7,, 

Here, prime denotes the derivative with respect to x = In a. Taking derivative of Qk = —k/(a 2 H 2 
with respect to x = In a, one gets 



ft' fc = -20fc(l + 



H 2 ) 



(29) 
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A. Noninteracting case 



Consider the FRW Universe filled with ECNADE and pressureless DM which evolves according 
to their conservation laws 



PA + 3il(l + ua)pa = 0, 



(30) 



p m + 3Hp m = 0, 



(31) 



where oja = Pa/ PA is the equation of state (EoS) parameter of the ECNADE model. 

Taking time derivative of the first Friedmann equation (|13p and using Eqs. (|23|) . (|'24p . ([25|) . 
(USD, d2TJ) and dM]), one can get 

H 1 



H 2 2 



ft fc -3(l-O0 

^/^A\3/2/ G c 8itG c aH 2 n A 
na V 7 n / V n G g 3n 4 7 n 



Substituting this into Eq. ([28|h one obtains 



na \ 7„ / 

Oa-1\/. G c 87rG c a# 2 A 



7n 



2 7 n 



a 



3n 4 



7n 



Putting Eq. (|32|) into ([29]) reduces to 

(l-O fe )-30 A 



na V 7„ / 

The deceleration parameter is given by 



J2_/^\3/2/ G c 8itG c aH 2 » A 

no V 7„ / V n G g 3n 4 7„ 



1 + 



-1 

# 2 / 



Replacing the term H/H 2 from (|32p into (|35p yields 



« = ^(l-^-3^) + -(^) 3/2 (27n-^ 



na V 7„ / 



G g 



8irG c aH 2 fi A 



(32) 



(33) 



(34) 



(35) 



3n 4 7 

Note that if we set a = (3 = then from Eq. fl26|) 7„ = G c /G g . Therefore Eqs. jM 
reduce to 



(36) 
and d36l) 



£2. = -J7 A 



0, + i<>A- D(3-^!(^ Vi 2 
na V G c 



(37) 
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2 na \G C , 



(38) 



which are same as those results obtained for the NADE model in HL cosmology (Jamil and Saridakis 
2010). 

Also for G c —> G g , Eqs. (|33p and (|36p recover the results obtained for the ECNADE model in 
the standard FRW cosmology (Karami and Sorouri 2010) as 



Q' A = fi A 



3(1 - n A ) - n k + 
n A -i 



7r, 



2 7 r, 



1 



na V 7„ / 
a# 2 ft A 



(39) 



g = -(l-fi fc -3fi A ) 



na V 7n / V 7n 3M 2 n 4 7n 



(40) 



The EoS parameter of the ECNADE model can be parameterized as (Huterer and Turner 1999, 
2001; Weller and Albrecht 2001) 



(41) 



Using ([30j) and flU}, the ECNADE density evolves as (Huterer and Turner 1999, 2001; Weller and 
Albrecht 2001; Copeland et al. 2006) 

PA 



PAo 



The Taylor expansion of the DE density around ao = 1 at the present time yields 

dlnp A 



In p A = In p Ao + 



din a 



I d 2 ln^A 
o + 2d(lna) 2 



(hi a) 2 + 



(42) 



(43) 



where the index denotes the value of a quantity at present. Using In a = — ln(l + z) ~ — z + 4r 



(44) 



for small redshifts, Eqs. (|42p and (|43p . respectively, reduce to 



In {p\/ PA ) 



In a 



-3(1 + w ) - 2^ lZ ' 



In (p A / PAo) dlnp A 



1 d 2 In p A 



o 2d(lna) 2 



2. 



Comparing Eqs. 



In a d In a 

and ([35]). one can obtain the parameters wo and wi as 

1 din p A 



w 



3 din a 



1, 



(45) 



(46) 
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1 d 2 In pa 

^1 = 77 







3d(lna) 2 

From Eq. (f3~T]) . the energy density of DM evolves as p m = p mo a ~ 3 - Now using (|24l) one can get 

-3 



(47) 



Pm n 
PA = 7^— "A 



Pm a 



(i - o fc - n A ) 



n A . 



(48) 



Substituting the above relation into ([46]) yields 



w 



ir^ A+ n A + fi' fc 



o A (i - n k - n A ) 



(49) 



Adding Eqs. {33} and ([Ml), we get 

ft A + J7' fc = (l-O fc -0 A ) 



O fe + 3^-A^ 3/2 



na V 7„ / 



x 27 r , 



G c MG c aH 2 Q A 



Go 



Substituting Eqs. ([331) and into (09]) yields 



w 



- 1 + ^„V7,„ 
8vrG c a^ 2 ft Ao 



3n 4 



^A W 2 



7n 



2 7n " 7^ 



(50) 



3n 4 7„ 

Note that in the absence of correction terms, i.e. a = (3 = and 7 n = G c /G g , Eq. ([ST]) yields 



(51) 



0J 



-1 + 



(52) 



3ra \G C . 

which is the EoS parameter of the NADE model in HL cosmology (Jamil and Saridakis 2010). 

Also in the limit of G c — > G g , Eq. (|5ip recovers the EoS parameter of the ECNADE model in 
the standard FRW cosmology (Karami and Sorouri 2010) as 

2 /O\ \i/2 



UJQ 



-1 + 



3n7n V 7„ 



x (27„ 



1 



Using Eqs. ([4"T]) and ([48]) . one can get 



+ 



3M 2 n 4 7no 



^A + «2 



(53) 



3Lft A ft 2 ' (i-O fc -0 A ) 



/ \2 
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Taking derivative of Eqs. (|33|) and (|34J) with respect to x = In a, one can obtain 



and 



o A = -(o A o' fe + o A o fc ) + 30 A (i -20 A ) 

rG r 2n 2 



+(o A -o A -i) 

40 A 



40 A + 16vrG c a 



+ (1-0/ 



Gg aH 3 r] 3 aHr] 3aH 3 r] 
G c 6n 2 \/H 



V afir? G g a#V / V fi 2 #77 / 



+ 



■rj_\ 



40 A 16vrG c a/3i? 5r) \ 
~~aHr~) ~ 3aH?>rf + 



(55) 



n' fc (i - 2n k ) - 3(o A o' fc + o A o fc ) + (n k - n' k ) 



G r 2n 2 



40 A + 16vrG c a 



G g aH 3 r] 3 aHr] 3aH 3 r/ 5 . 
y40 A G c 6n 2 wil 77 \ 
XaHr] ~ G^alPrf) KH 2 + H~r~]) 
16vrG n a /3H 



-0 4 



40' 



+ 



_5rj_\ 



aHr] 3aH 3 r] 5 \H 2 Hr]J. 
Using Eq. (j5T]l . the above expressions for A and 0'^ can be rewritten as 

A = — (OaO'j, + A Ofc) — 30a(0a — A — 1)ujq 

+(1 - o A ) h(n A + A ) + a] , 



(56) 



(57) 



where 



= 0' fc (l-20 fc )-30 A (O fc -0' fc H 



.4 



Using Eqs. ([35]) and r) 



-Oi 



3(0 A + 0' A ) + A 



G c 6n 2 



40', 



/ 40 A 

\aHr] G g aH 3 r] 3 J \H 2 HrjJ aHr] 

167rG c a (3H 5r) \ 
~3aH 3 r] 5 \lP + fir//' 

z l/a, one can rewrite f|59[) as 

A = -90 A (l + g)(l+o;o) 

+ ^(^) 1/2 fl7 + 4(O fc + 2,) 

+3(7-40aVo-^(^) 1/2 

32irG c aH 2 /Oa\ 5 / 2 ' 
3an 5 OA v 7n ' 



(58) 



(59) 



(60) 
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Adding Eqs. ((37} and ([35]). we get 



-(n A n' k + n' A n k ) + n' k (i-2n k ) 
-3(n fc - nj. + o A - ^ A - 1)^AW 
+(i - n A - n fc ) k^A + o A ) + a 



Substituting Eqs. ((33 ]) . ((5Djl . ((37]) and (|6Tj l into Q5JJ) yields 

wi = (1 + w ) 1 - fife - 3(1 - Ao Vq 
which can be rewritten by the help of Eq. (|60|) as 



+ 



An 



30 Ao 



vl+W ) 



1 - O feo - 3(1 - ^a )o;o 



1 /^An\ 1 / 2 

3n V 7„ Q / 
17 + 4(fi fco + 2q ) + 3(7 - 4fi A( >o 
8 /Q Ao y/2 327rG c qg 2 / 17 Ao \5/2 
n V 7„ / 3n 5 Ao V 7 no 



(61) 



(62) 



(63) 



B. Interacting case 

Here, we consider a case in which the ECNADE and DM interact with each other. This causes 
the energy conservation law for each dark component not to be held separately, i.e. 



PA + 3H(1 + uj\)pa = -Q, 



(64) 



Pm + 3Hp m = Q, 



(65) 



where Q = 3b 2 HpA (Pavon and Zimdahl 2005) stands for the interaction term with coupling 
constant b 2 . For Q > 0, there is an energy transfer from ECNADE to DM. The recent observational 
evidence provided by the galaxy cluster Abell A586 supports the interaction between DE and DM 
(Bertolami et al. 2007, 2009; Abdalla et al. 2009). 

Taking time derivative of the first Friedmann Eq. flT3]) and using (p3"1) . (giD , (p31) . ((261) . (pTj) and 
(1631) yields 



H 
W 



n k - 3(1 - n A ) + zb 2 n A 



X ( 2 7n - g- 



J_fttA_y/ 2 

na V 7„ / 



87rG c aH 2 Cl A 
3n 4 7n 



(66) 
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Substituting this into Eqs. ([28)) and (]29j) one can get 



3(1 - n A ) - 3b 2 n A -n k + 



na \ 7 n / 



'^A-l\/ G c 8nG c aH 2 n A 
x I I 27„ - — - 



7n 



(So 



3n 4 7„ 



and 



J7i — r^t 



(i - n k ) - 3n A - 3b 2 n A + —(^) 3/2 



na \ 7 n / 



/ G c 87rG c aH 2 n A 

x 27 n - —4 

V G g 3n 4 7 n 

Replacing the term H/H 2 from (|66p into (|35p yields 



From 



2 na V 7„ / V 



lnG c aH 2 Q A 
3ra 4 7„ 



— 6 Oa- 



(67) 



(68) 



(69) 



one can obtain 



PA 



'1 — — J7a) 



dln/9 A _ p' m ^ A 



din a p m Q m fl A 
From Eq. (|64p and using ()41|) . the interacting ECNADE density evolves as 

P±_ _ Q -3(l+o;,)-a;i+6 2 ) e 3L l ;iz 
PAo 

2 

Using In a = — ln(l + z) ~ — 2 + 4r for small redshifts, Eq. ([72]) reduces to 



In (pa/Mo) 



-3(1 + o;o + 6 2 ) — -u)\z. 



(70) 



(71) 



(72) 



(73) 



In a ' ' " ' " ' 2" 

Comparing Eqs. (145p and (j73|) . one can obtain the parameters wo an d u)\ for the interacting case 
as 



too 



1 din pa 
3 din a 



(74) 



Id 2 In 



^1 = 77 



PA 



3d(lna) 



(75) 
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Substituting Eq. (|7T|) into ([71]) and using ([65]) one can get 



u; 



0' A 



+ 



n A (i-n k -n A ) 



1 — Q k ~ &-a' o 



Using Eqs. ([70]) and ([75]) one can obtain 



3b 2 n' A 



l-Oi 



ft A (l-n fc -n A )2 + A 



/ \2 



A ' (l-Ofe-fi A ) ' (1 



Adding Eqs. ([67]) and ([68]) gives 



n' A + n' k = (i-n k -n A ) 



n k + sn A -^m 3/2 

na V 7„ / 



x 27„ 



G c MG c aH 2 n A 



Go 



3n 4 



7n 



-3b 2 n A (n k + n A ). 



(76) 



(77) 



(78) 



Substituting Eqs. ([67]) and ([78]) into ([76]) yields 

2 /^Ao\V2 



(j 



-1 + 



3rJ7„ V 7ri0 
/ Gc ^G c aH 2 n Ao 



Go 



3n 4 7 no 



(79) 



which in the absence of interaction, i.e. b 2 = 0, reduces to Eq. (|5ip . The EoS parameter ujq 
versus the interacting coupling parameter b 2 is plotted in Fig. [1] Figure shows that in the absence 
of interaction (b 2 = 0), ojq = —0.85 and cannot cross the phantom divide line. However, in the 
presence of interaction, b 2 ^ 0, the phantom EoS ljq < — 1 can be obtained when b 2 > 0.15 for the 
coupling between ECNADE and DM. This value for coupling constant b 2 is consistent with recent 
observations (Wang et al. 2005a, 2005b). Also the phantom divide crossing is compatible with the 
observation (Komatsu et al. 2011). 

Taking derivative of Eqs. ([67]) and ([68]) with respect to x = In a and using ([79]) . one can get 



-(n A n' k + n' A n k ) 
sn A (n A - n' A - i)(gj + b 2 
+(i - o A ) h(n A + n' A ) + a 



-6b 2 n' A n A , 



(80) 
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and 



«2 



n' k (i - 2n k ) - 3n A (n k - n' k )(cj + b 2 ) 



a;> 



(81) 



where ^4 is given by Eq. ([59]) . Using Eqs. ([25]) . ([35]) . (J66J) , ([79]) and 77 = 1/a, one can rewrite ([59]) 
for the interacting case as 

A = -9O A (l + g)(l+c + & 2 ) + ^(^) 1/2 



17 + 4(ft fc + 2q) + 3(7 - 4fi A )wo + 216 2 
_8_/^A\i/2 32irG c aH 2 /OaW^ 
na V 7„ / 3an 5 f2 A Wn ' 



Adding Eqs. ([SOB and JST]) gives 



-(n A n / fc + i2^n fc ) + f2' fc (i-2n fc ) 

-3fi A (wo + & 2 )(^fc - Ofc + A - Q A 



1) 



+(1 - Oa-^a 



3(ft A + A ) + A 



-36 2 (^a + ^k^ A ) - 6b 2 n' A n A . 



Replacing Eqs. <[S7 j> . ([75]) . ([50]) and (|H5| I into (J77J) gives 

o 2 ft Ao 



(1 - ^Ao - &k ) L 

+6Q Ao u} - 2£l ko 
+(30a ^o - tt ko ) 

b 2 n Ao 

(1 - ^A - ^fc 
6 2 ^A 



3(w + Zf 



L 3 



+ 



(1 - J7ao - flk ) A ° ° + 3 



x 3w (l -^Ao) + 3o 2 +O feo 
+(w + & 2 ) 3(w + 6 2 )(20 Ao - 1 
+1 - 2Q ko - 6b 2 n Ao 



+ 777^ + 1, 



3^A 



(82) 



(83) 



(84) 



which in the absence of interaction, i.e. b 2 = 0, yields Eq. (|63|) . The EoS parameter u\ versus 
the interacting coupling parameter b 2 is plotted in Fig. [2j Figure shows that in the absence of 
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interaction (b 2 = 0), oj\ = 0.07 and the EoS parameter ui\ increases when the interacting coupling 
parameter b 2 increases. This is in agreement with the result obtained by Jamil and Saridakis 
(2010). 

IV. CONCLUSIONS 

Here, we investigated the ECNADE scenario in a FRW Universe in the framework of HL gravity. 
We considered an arbitrary spatial local curvature for the background geometry and allowed for an 
interaction between the ECNADE and DM. In both the regular and interacting case we obtained 
the deceleration parameter as well as the differential equations which determine the evolution of 
the ECNADE density parameter. Finally, using a low redshift expansion of the EoS parameter of 

DE as uj\(z) = ujq + oj\z, we calculated ojo and oj\ as functions of the DE and curvature density 

G 

parameters, VL\ {) and ttk respectively, of the running parameter A = (2^- + l)/3 of HL gravity, of 
the parameter n of ECNADE, of the interaction coupling b 2 , and of the coefficients of correction 
terms a and (5. 

Interestingly enough we found that in the presence of interaction between ECNADE and 
DM, the EoS parameter ojq of ECNADE in HL gravity, can cross the phantom divide line which 
is compatible with the recent observations. The interaction between ECNADE and DM can 
be detected in the formation of large scale structures. It was suggested that the dynamical 
equilibrium of collapsed structures such as galaxy clusters would be modified due to the coupling 
between DE and DM (Bertolami et al. 2007, 2009; Abdalla et al. 2009). The idea is that the 
virial theorem is modified by the energy exchange between the dark sectors leading to a bias in 
the estimation of the virial masses of clusters when the usual virial conditions are employed. This 
provides a probe in the near Universe of the dark coupling. The other observational signatures on 
the dark sectors' mutual interaction can be found in the probes of the cosmic expansion history by 
using the type la supernovae (SNela), baryonic acoustic oscillation (BAO) and cosmic microwave 
background (CMB) shift data (Guo and Ohta 2007; Feng et al. 2008; He et al. 2009; Honorez et 
al. 2010). 
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FIG. 1: The EoS parameter ljq of the ECNADE in HL gravity, Eq. (|79|) . versus the interacting coupling 
parameter b 2 . Auxiliary parameters arc: n = 2.716 (Wei and Cai 2008b), a = —7.5, 7„ = 15 (Karami and 
Abdolmaleki 2010c), Q Ao = 0.728, fl ko = -0.013 (Komatsu et al. 2011), A = 1.02 (Dutta and Saridakis 
2010), G c /G g = 2/(3A - 1), H = 74.2 Km S" 1 Mpc" 1 (Riess et al. 2009) and M p 2 = 8ttG s = 1. 
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FIG. 2: The EoS parameter u>i of the ECNADE in HL gravity, Eq. (|84|) . versus the interacting coupling 
parameter b 2 . Auxiliary parameters as in Fig. [TJ 
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